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A Perturbation Method for Studying Thermal Convection Problems 


By 


Bin Oxat 


Geophysical Institute, Faculty of Science, Tokyo University, Tokyo 


Abstract 


When a layer of fluid is heated uniformly from below, and the Rayleigh number 
exceeds a critical value, convection takes place in it in a regular cellular pattern. A per- 
turbation method is presented here for determining mathematically the pattern and 
amplitude of this steady convective motion. The essential point of the method is to 
expand functions which describe the velocity and temperature field in the fluid in a power 
series of a parameter ¢, while the Rayleigh number is put as a product of the critical 
value times (1--e?). A set of inhomogeneous equations thus obtained can be solved by 
the perturbation method which is in use in non-linear oscillation problems. In the two- 
dimensional case, the slope of heat transport curve becomes steepened abruptly when the 
Rayleigh numher exceeds a critical value. Another problem which can be dealt with by 
this method is that of convection within a sphere. This forms an extension of CHANDRA- 
SEKHAR’s linearized stability theory. 

Furthermore, a study is made of the steady thermal convection in a two-dimensional 
fluid layer when it is heated uniformly from below under a simultaneous constraint of non- 
uniform temperature given on its upper surface. Mathematically this is an application of 
the method mentioned above to a problem with inhomogeneous boundary conditions. It 
_was found that the site of spontaneous convection cells is governed by the surface tem- 
perature disturbance having a critical wave length. Surface temperature disturbances 
having much longer or shorter wave lengths play very little part in this, while those 
having wave lengths close to the critical one have influential effect in determining the 


general feature of fluid motion. 


§1. Although the earth’s mantle behaves as 
a solid for forces of short duration such as 
those occurring in earthquake waves and in 
earth tides, it may well behave as a fluid for 
those of very long duration, such as suggested 
by approximate isostatic adjustments of large 
topographic features. If this is actually so, 
there is a possibility that the mantle may be 
subject to a very slow convection (PEKERIS, 
1935; Hates, 1936; Buniarp, 1950). Of course, 
some difficulties arise in supposing convection 
currents in the mantle. Brrcw (1951) has 
pointed out that the currents, if they take 
place at all, can hardly extend up to the 20° 
discontinuity. Two recent discoveries appear, 
however, to support the hypothesis of convec- 
tion in the mantle. One is concerned with 
the heat flow from within the earth. The 
measurements which have been made up to 
now in this connection (BuLLARD, 1954; BuL- 


LARD, MAxwELL and RevELLE, 1956; Herzen, 
1959) appear to indicate that the heat flow in 
oceanic bottoms does not differ significantly 
from that in continents. The source of heat 
is reasonably supposed to be due to radio- 
activity within the crust, which is 35 km thick 
on the continent and 5km under the ocean. 
Where does the oceanic heat come from? If 
radioactivity is distributed through too great 
a depth, the temperature near the bottom of 
the radioactive layer will rise above its melt- 
ing point, and this is inconsistent with the fact 
that S waves can be propagated at these 
depths. This difficulty may be avoided if the 
transport of heat by convection in the mantle 
is assumed. The other discovery is concerned 
with palaeomagnetism. Studies of palaeomag- 
netism (BLACKETT, 1956; Runcorn, 1955; 
Creer, IRvING and Runcorn, 1954) have re- 
juvenated interest about the classical hypo- 
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thesis of continental drifts. But some reason- 
able forces which drive continental blocks in 
motion must be looked for if the hypothesis 
is to become more than a mere hypothesis. 
Here again, these forces may be assumed to 
be due to convection currents within the 
mantle. 

Although there are evidences like these in 
favour of convection currents in the mantle, 
nothing definite can be said at present about 
the existence of such large scale motions it- 
self. This is partly because we do not have 
enough knowledge about the earth and partly 
because mathematical studies of convection so 
far developed are far from being satisfactory 
for the purpose. There are many funda- 
mental problems of thermal convection to be 
solved. How much heat is transported after 
convection once takes place in a layer of fluid 
which so far has been at rest? What would 
the distribution of temperature become in a 
convective layer? In 1958, MaLkus and VsERonIs 
presented a perturbation method by which the 
form and amplitude of cellular convection 
can be determined. Their intention was to 
complete RAYLEIGH’s first-order solution of the 
problem. At about the same time the present 
writer found a similar method independently. 
It was when he had finished his calculations 
up to §8 of the present article, that he was 
informed of Ma.kus and Veronis’ paper by 
Dr. CHANDRASEKHAR Of Chicago University (Oct. 
24, 1958). A part of the results presented 
here is similar to those of Markus and 
Veronis. In the following, stress will be laid 
upon the writer’s own method of analysis 
which is more or less different from theirs. 

In §2, the fundamental equations and bounda- 
ry conditions of the problem are stated. In 
§3, the perturbation method of solutions for 
the finite amplitude fields of temperature and 
velocity is described. The method can be 
applied to other problems, of which a simple 
example is given. In §4, this method is ap- 
plied to a two-dimensional convection problem 
(roll case). In §5, the relation between the 
heat transport and the mean temperature 
gradient is studied. In §6, the fundamental 
equations and boundary conditions are describ- 
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ed for the convection within a sphere and in 
§7, analysis is carried out on this problem. In 
§8, convective fluid motions which take place 
within a slowly rotating sphere are considered. 
Up to this point the boundary conditions are 
homogeneous. Hereafter inhomogeneous 
boundary conditions are dealt with. The 
reason why the inhomogeneity is brought up 
is as follows. When the temperature dif- 
ference between the upper and lower surfaces 
exceeds a certain critical value, convection 
takes place in a regular cellular pattern, but 
theory does not distinguish whether horizontal 
pattern of convection cells is rectangular, 
triangular or hexagonal. Thus we cannot 
answer a question such as; Why does one of 
these patterns take place preferentially, when 
more than one solution is possible at the point 
of instability? Even when the motion is taken 
to be two-dimensional, one cannot predict 
where the site of upwelling will be located in 
the fluid. The motion could be upward any- 
where, so long as the fluid is uniform hori- 
zontally. But such a perfect uniformity 
cannot be imagined in nature. It is of im- 
portance therefore to know the influence of 
inhomogeneity. Besides, in the case of con- 
vection within the earth’s mantle, considera- 
tion must be paid to the earth’s crust of 
non-uniform thickness. §9 and §10 correspond 
to §2 and §4 for homogeneous case respec- 
tively. In §11 the results obtained are given. 
Finally it is admitted that there are many 
related geophysical problems out of the scope 
of the new method presented here. This is 
merely one step toward generalisation of the 
classical linear stability theory. 
§2. First of all, a two-dimensional steady 
convective motion in a fluid as shown sche- 
matically in Fig. 1 will be dealt with. The 
fluid is supposed to be enclosed between two 
horizontal conducting surfaces, which are h 
apart. The lower surface is kept at a tem- 
perature 7) and the upper surface at a lower 
temperature (Z»-8h), the mean temperature 
gradient being -8, (8>0). The fundamental 
equations are the equation of continuity (1), 
the equation of motions of Navrer-Stoxgs (2), 
(3), the equation of thermal conduction (4) 
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x (8) 


F: = Te 


Fluid layer heated from below. 


Bice l: 


and the equation of thermal expansion (5): 


a Lae (1) 
aioe Ot) 

o( uf ay ee +uPr7u (2) 
Cie es Nee Fy ert ary (3) 
u or +05 ser (4) 


0=po.{1—a(T—T»)}, (5) 
where uw and v are the components of fluid 
velocity in the direction of x and y respec- 
tively, J is the temperature, p is the pres- 
sure, g is the magnitude of the gravitational 
acceleration, yu is the viscosity of the fluid, « 
is the thermometric conductivity, o is the 
density, and a the linear coefficient of thermal 


expansion. Using the stream function ¢, u 
and v may be rewritten as u= Hs and 
0g 
v=—— 
Ox 


At this point it is of value to non-dimen- 
sionalize the physical quantities appearing in 
the equations. This will simplify subsequent 
mathematical manipulations and help us gett- 
ing pertinent physical parameters for the 
problem. 

Let T=Bh(—y+1), w=hé, y=hy, o=Ky, 
(6 ) 
and after these transformations, all primes 
will be dropped out. Then (2) and (3) are 
cross-differentiated and the difference taken 
in order to have pressure eliminated. The 


resulting relations between ~ and r are 


LGB 8-2 (8228 


On? §=0F0x 
= Or 4 
aR eee (7) 
Of Oc Ob Or Ow A 
Oy (OF. O6c dpe idee. ohee 2 
4 
Ra see is a non-dimensional parameter 


which is called the RAYLEIGH number and 


es is the PranpTtL number. So far ail the 


non-linear terms have been retained, but, 
hereafter, for the sake of simplicity, the 
inertia terms will be ignored, as it is those 
non-linear terms appearing in (8) that are 
more essential in the problem. Also, in the 
case of convection within the earth’s mantle, 
the inertia terms are negligible because the 
pertaking viscosity is enormously high. Thus 
the final equations become 


Op Or Of Or Op _ 
SF Ot On OF ae ) 
Py Le (10) 
Ou 


The conditions at the boundaries (y=1, y=0) 
are 
Oh Ob 
Oy? Ox? 
following RayteicH. This means that tem- 
peratures are kept constant, the vertical 
component of velocity vanishes and stress 
does not exist at the boundaries. In other 
words, the boundaries are free. More realistic 
are of course the conditions of ‘rigid-rigid’ or 
‘rigid-free’, but these are difficult to be taken 
into account. Essential points of the problem 
are not lost by adopting the above-mentioned 
unrealistic conditions. 


=(0, r=0, (11) 


Oe 


§3. A convective motion takes place in the 
fluid layer when the mean temperature gra- 
dient 8 (or R) reaches a critical value fo (or 
Ro). A further increase of mean temperature 
gradient produces more intense field of motion. 
Therefore functions # and tr describing the 
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field must appear suddenly when § reaches 
8). In other words, # and r are not analytic 
at B= (Fig. 2). Consequently one cannot 
expect solutions in the form 


B=B6,(1+0) | 
b=0-dit6?-d2+0'-b3s+-:- i 
t=60-71 +0? -t2+03-t34+:-: 


which are all analytic. If ~ and t could be 
written in the form mentioned above, # and 
t could be extended to include the state 
having the temperature gradient less than the 
critical one, and this will mean that a convec- 
tive motion could occur even before the mean 
temperature gradient § reaches the critical 
value So. This is absurd and does not re- 
present actual cases. Therefore another form 
of solutions has to be looked for. The weakest 
point of the expression B=§,(1+6) in (12) is 
that 8 can easily become less than $0, while 
physically 6 must always be positive. This 
can be avoided if 6 is put to be equal to &. 
But if all the 6’s in (12) are made equal to 
&, the situation will still be the same, for 
there is no reason to suppose that € does not 


Wor G 


(12) 


A, P 


Fig. 2. Relation between mean temp. grad. and 
the function describing the field. ¢, (c) is not 
analytic at the point where B=). 


become imaginary and the absurdity remains. 
The difficulties will disappear if one could find 
solutions in the following form 

(13) 


B=8(1+ &) 

P=Edi terdet+E %p3+:-- 
T=EtT1 + E*t2+6%r34+-:- 

Of course, here again 8 becomes less than 

if € is taken as pure imaginary. But ~® and 

7, then, become complex functions which have 
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no physical counterpart in reality. This means 
that » and r as given in (13) are good enough 
for describing the field. It is now required 
that expressions (13) satisfy the equations of 
motion, (9) and (10), for all values of € less 
than some maximum &m. The coefficients of 
each power of & obtained on inserting the 
expressions (13) into the fundamental equa- 
tions (9), (10) must vanish individually and 
the resulting series for / and tr must converge 
if relations (13) are to represent a satisfactory 
solution to the problem. 

The sequence of linear inhomogeneous equa- 
tions in powers of & obtained on inserting 
relations (13) into equations (9) and (10) is 


At €'d4=Ro(1+€?) o Dene (14) 


Roe i 0x i 5 Me i rae i 
Of te oe oe 


+ Semper. (15) 
ct i 

Using # andr, the original quantities uw, v and 

T are expressed as follows: 

Ob 


== 
Oy 


aoa T=hpi PE) Hey 
(16) 


v=— 


iS 
h 


>| 


where 2 and y are non-dimensional. 

The expansion method described here is 
also applicable to other problems, in which 
dependent variables are single-valued and the 
domain of these variables is bounded. As a 
simple example, the problem of one-dimen- 
sional simple harmonic oscillator is dealt with. 
The governing equation with regard to the 
velocity v is 

» + w?2=0 ) 
which, of course, can be integrated easily. 
But here, the above-mentioned method of ex- 
pansion is followed. In (17), x is upper- 
bounded and has some maximum value, a 
(amplitude of oscillator), as shown in Fig. 3. 
The situation near x=a (or e=—a) in Fig. 3 
is quite the same as the neighbourhood of 


B=) in Fig. 2. Therefore v and x are put 
as 


e=a—-&, V=AiE+ A624 Az63+.--, 


(17) 


(18) 
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By changing the independent variable in (17) 
from 2 to &, (17) becomes 


ay 


Fig. 3. Relation between velocity and displace- 
ment of a linear oscillator. v is not analytic 
where x«=+a. 


v dv ee ae 
26 d& +a?(a—é&2)=0. 


—_— (19) 
Explicit determination of Ai’s need not be 
mentioned, as it is quite elementary. Finally 
it is shown that the solution thus obtained is 
wo 


no other than SU+ Stak. 


§4. In this section, the sequence of equa- 
tions, (14), (15) will be solved. 
[€]: The first-order equations are 


071 


Ath; =R 
di =Row 


°e 


(20) 


* pare (21) 
These are concerned with the stability theory 
and have been exactly solved in terms of 
_trigonometric functions (Rayueicu, 1916). A 
complete solution is 


¢1= A: sin zy sin ame (22) 
i= oS sin zy cos a (23) 
Ro= “txt (24) 


This is obtained by assuming a solution in 
the form e’”’“*sin mzy and then making the 
value of R (or the mean temperature gradient) 
as small as possible. Thus Ro=Min R(m, 1), 


m, n \integer) 
where n is an integer. An explicit value of 
A; in (22) and (23) will be obtained in the 


third-order stage [€*]. Finally it is noted at 
this point that the origin of the coordinates 
is so chosen as to make w=0 at «=0. 


[e*]: Pa Roe* (25) 
Sis Ors _ 00s Ori Os ore, (06) 


Gy Oz Ox Oy | Ox 
Inserting (22) and (23) into (25) and (26) and 
solving the equations thus obtained, we get 


the following solution to the second-order 
equations: 


do= As sin zy sin 5 (27) 
T= V2 A; sin zy cos ——— ge sin 2zy. 
3x Va 247 
(28) 


Here again, the solution contains a new 
arbitrary constant A,, which is determined 
explicitly in a later stage. 


LG2le rds=R 08s ORE Or. (29) 
°' Ox Oa: 
Od, OT. Ode Or1 — Opn Or Ode Orr 
Oy Ox Oy Ox Oz Oy Ox Oy 
4a ae = =*rs (30) 
we 


Inserting (22), (23), (27) and (28) into (29) and 
(30), we get the following equation: 


2 2 
Pn Roo = an Ae —1 sin mysin a 
SR Sain (31) 
ei sin 3zy sin a5 x. } 


On the right hand side of (31) there are two 
inhomogeneous terms. The one which has 


the form of sinzysin 73 will generate a 


particular integral which does not satisfy the 
boundary conditions. Therefore we have to 
choose the value of Ai so as to eliminate 
the resonant inhomogeneous term in (31), (a 
method used in non-linear-oscillation problems): 


Rox* Lert \=0 


Ai=V/24, —v/24, 0 (33) 
Thus the determination of A: has been made 
which was impossible in the earlier stage of 


(32) 
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[é]. Of the three values in (33), —V/ 24 is 


discarded as it coincides with 7/ 24 by another 
choice of coordinates. Also A:i=0 must be 
discarded. This represents no other than the 
state of unstable rest. 


Rye Ot 


[é4]; d= +R, “On (34) 
Od Ors | Ops Ore oe Orr _ Odi Ors 
Oy Ox Oy Ox Oy Oz Ox Oy 
_ Ope Ove _ Obs Or ee 2 
Ox Oy Ox Oy ° Ox ae) 


Rearranging the equations into a single one 


[é] a ws 
$1 = 4.898 sin zy sin ve ze 


= — 0.735 sin zy cos 


= 


V 2 
[E*] (¢ ¢2=0 
= —0.318 sin 2zy 


5 


t3=0.675 sin zy cos 


V2 
[é*] 


—*_4—0.116sin 37zy COS—: 
V 
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and eilminating the resonant inhomogeneous 
term, we get the following relation as in case 


of [€°]: 
Ro Az 2 Ay )=0 36 
9 7 (2 3 (36) 
A2=0 (37) 


From A.=0, it is expected that the solution 
in the form (12) will vanish identically, for 
(12) is no other than (13) having the relations: 
di=ds=---=u=ta=---=0. 

Explanation of the following procedures will 
not be necessary. The results are just written 
down up to the stage of [6&7]: 


—a+0.019 sin 3zy sin a 


VA 


re 


¢4=—0.058 sin 2zy sin Y 2 xx+0.001 sin 4zy sin Y 2 zz 


t4=0.320 sin 2xy—0.025 sin 4zy+0.070 sin 2zy cos Y 2 xx 


le | 3=0.385 sin zy sin 


—0.011 sin 4zy cos Y 2 xx 


[€5] os=—0.130 sin zy sin 


vA 3 
—0.001 sin 3zy sin ue 2 Tx 


Ts = —0.666 sin zy cos i = 


+0.004 sin 3zy cos a8 2  — 


[é°] 


e=—0.323 sin 2zy+0.034 sin 4zy—0.001 sin 6zy—0.1500.Ae sin zy cos ———= 


-7-+0.070 sin zy sin 3v2 za —Q.006 sin 3zy sin ——= 


>= 5%—0.047 sin xy cossV. 2 


do= Ao sin zy sin ae «+0.041 sin 2zy sin Y 2 mx—0.014 sin 2zy sin 2 2 nz 
—0.001 sin 4zy sin Y 2 xa 


Z Jr? 


zet-0.157 Sin 37y ‘cos —— 


2 Vr" 


0.009 sin 5zy cos ——— 


73" 


yr 


—0.119 sin 2zy cos Y 2ze+0.005 sin 2za cos2 Y 22x 


+0.017 sin 4zy cos Y 2 xz+0.001 sin 4zy cos 2 2 xx 


[€"] 
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,;= A; sin zy sin 


Vr 


7=—0.130As6 sin 2zy+-(0.666—0.1500.4;) sin zy cos 


+0.061 sin zy cos 


—0.176 sin 3zy cos - 


+0.013 sin 5zy cos rs «+0.001 sin 5zy cos 


—0.020 sin zy sin 


V2 


3y —zx—(0.001 sin zy cos 


2 


V2 


x—().001 sin 3zy cos 


zx+0.003 sin 3zy sin 


BY 2 
2 
va Tes 
oo 
3Vv z TX 
Es 
an 2 a 


Vitg 
ee 


(38) 


Explicit determination of Ae and A; has not yet been made, though As is expected to vanish. 
(38) is rewritten to show the coefficient of each trigonometric function in a power series of &: 


yp: 


sin zy sin—“—2: 4.8986 +0.38563— 


v2 


Sawai ee. _ 0,070°—0.0208 


V2 


sin 2zy sin VY 2 xz: 


sin 2zy sin Len 2: 


sin 37zy sin fy: 


s 2 
sin 3zy sin a. 


yo o2 


sin 4zy sin VY 2722: 


— 0.0584 +0.041€° 


—0.014é° 


0.019€%—0.006€° +-0.003€" 


—0.001€° 


0.001é*—0.005€° 


sin 2zy: —0.318€?+0.320E64—0.3236°— 
sin 4zy: —0.025&4+0.034E* 


sin 6zy: —0.001E° 


Sil COS = — 


75? 


ea ae of <= 
sin zy te 


57 
sin zx COS V Pag Mae 


sin 2xy cos VY 2 za: 


sin 2zy cos 2V/ 2 a: 


sin 3zy cos 


7 a 


in 3 os eset 
sin 37y C Jo : 


sin 47ry cos VY 2 rw: 


—0.735€ + 0.6753 — 


0.666€° 


+(0.666—0.1500.A7)é7 


—0.047&°+0.061E" 


—0.001€7 


0.07064—0.119€* 
0.005€° 


— 0.116? +.0.1576°— 


0.004€°—0.001E" 


—0,011€*+0.017€° 


0.1766" 


0.130. A%E" 


0.1306°+ A,é*+ A,é" 


—0.1500A,° 
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sin 4zy cos2V 22x: 0.00164 


7 
j a — .009E5—0.013E" 
sin 57 cos e G 0 
one COn ee OLE! 
sin 5zy rap ; 


§5. The heat transport increased by convec- 
tion is easily calculated from the temperature 
gradient at the boundaries: 


dT ( a.) 
=— =a = (ou = 
Heat flow Gi FE kB)(1+ €?) Fee 
=hkB(1+6?)(1+1.996 -1.70€2+1.6163—- - -) 


(39) 
where the bars above TZ and ct denote an 


average over a horizontal line. Fig. 4 shows 

Heat flow/pp (1) 
25 1 
0 


=| -0.5 0 0.5 
Temp. grad, in excess 
of the critical value 


Fig. 4. Heat transport vs. temp. grad.. 
Dashed line represents the relation due 
to conduction alone. 


§6. Unless the fluid motion is restricted to 
two-dimensional, an infinite number of steady 
finite amplitude solutions (having different 
horizontal plan-forms) are obtained which 
formally satisfy the equations. 

Another case which also gives a unique 
solution as in the roll case is convection 


ae eS ee ee eee 


st Ta 
= MM 


the relation between the heat transport and 
the mean temperature gradient. As an 
illustration of convective patterns, isothermals 
and velocity vectors of the fluid are shown in 
Fig. 5, when the mean temperature gradient 
is 10% in excess of the critical value. The 
whole feature changes a little when the mean 
temperature gradient is 40% in excess of the 
critical value (Fig. 6). 


Fig. 5. Mean temp. grad. is 10% in excess of 
the critical value. 
(1) Heat transported from above. 
ordinate is the value due to conduction alone. 
Dashed line is the averaged value. 


Unit of the 


(2) ; Temperature and velocity distribution 
within a half of the cell. Ts and Tp, are 
temperatures of the surface and the bottom 
respectively. 


(3) Heat transported from below. 
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(1) 


el IT 


>> => 
—2Bu-r+u-gradt=xp?r . (42) 


In (42), the non-linear terms have been 
retained. Other equations governing 
the field are 


div u=0, (43) 


0=—grad p—o gradV+ ovp*u, (44) 


where the external field gradV is that 
due to the action of gravity in a 
homogeneous sphere, 
4 se : 
grad V=yrp0Gr. (45) 
G denotes the constant of gravitation 
and #o is the mean density. 


By replacing p by p=o(1—aT), we get 
from (44) 


Fig. 6. Mean temp. grad. is 40% in excess of 
the critical value (cf. Fig. 5). 


within a sphere, The following is an exten- 
sion of CHANDRASEKHAR’S linearized stability 
theory (1952). 

-Following CHANDRASEKHAR, the sphere is taken 
to be a homogeneous one of radius R in equi- 
librium under its own gravitation and with a 
uniform distribution of heat sources such that 
in the absence of conduction the temperature 
at each point will rise at arateof Q. Inthe 
steady state, the temperature distribution, 


T(r), inside the sphere will be governed by 

cp?T=—Q (40) 
The solution of (40) appropriate to the prob- 
lem on hand is 


a ‘ (41) 


T Har), B= 


In writing the solution in the form (41), we 
have assumed that T=0 at r=R; this entails 
no loss of generality. Let the temperature 
distribution in the convective state be given 


by T=T+r. The equation governing 7 is 


> > 
0=—grad ©+ytr++yp2u (46) 
1 
where hae nes V— — By (2R2r2—r4)} 
0 4 
(47) 
= 4 0G 
T= SP Qa, J 
The physical quantities appearing in the 


above equations are non-dimensionalized as in 
§ 2: 


a = 


CPR u= mur =RP. (48) 


After these transformations all primes will be 
dropped out. The following are the funda- 
mental equations in this problem: 


div u=0 (49) 
ru +u-grad t=f'r (50) 
0=—grado+-Goretpia ©”) 
6 
where ess : 
KY 


C corresponds to the RayLeiGH number R in 
§2. After eliminating ® from (51), we get 


r(ru)+ <p =0 (52) 


10 Bin 


h 1e=| 97% a, 2 
where = ety a? Ap ley 
Le tO ae 0 1 st 
~ sin? 00 at 00 i sin2d 0¢? 
and u=(u, v, w,) are the components of velo- 
city in the spherical polar coordinates. 
On the spherical boundary surface (r=1), it 
is required 


(53) 


1) that the surface temperature is fixed: 


T=0 (54) 
which, using (52), becomes V@w)=0, ©5) 
2) that w vanishes identically: FU=0, D6) 


and 
3) that, in case of a rigid surface, from the 
equation of continuity 


0 

wn =0, 57 
ar Ow (57) 
or, in case of a free surface, also from the 
equation of continuity 


2 


ne (ru)=0 . (58) 


But either of these conditions (57) and (58) 
makes following calculations extremely dif- 
ficult. Even in case of the linear stability 
problem, tedious procedures of variational 
method are necessary. Therefore, for the 
sake of simplicity, we impose the following 
artificial conditions: 


pone Mt my 9 

Pia ol, eae = ) = (00) 
These are reduced respectively to 

(4-+2)=0 (61) 

( o +=)w=0, (62) 


Using the equation of continuity, (61) and (62) 
become 

V(ru)=0. (63) 
From (59) and (60) it is seen that (63) is a 
condition between ‘rigid’ and ‘free’, which is 
also verified from (72) in §7. 


To sum up, the boundary conditions are (55), 
(56) and (63), 


OKAI 


§7. The perturbation method used in this 
section is quite similar to that used in § 3 and 
the solution is expected to be expanded in the 
following form 


C=Cil+é) 


w= Eur + E+ EUs: (64) 


tT=Et14+ 6t2+ 67734 --- 


where C> is the critical value of C and cor- 
responds to the critical intensity Qo of heat 
sources. 

Inserting (64) into (50) and (52) we geta 
sequence of linear inhomogeneous equations 
in powers of €. The following procedures 
are quite the same as in §4. 

{€]: The first-order equation is 

V°(ru1)=CoL*(ru). (65) 
Now we may assume for the solution of (65) 
the following form 

r= Wr) Y, ¢), (66) 
where Y is the spherical harmonic of the 
first order, as the lowest order gives C the 
least value. The equation for Wi(r) becomes 


2 oe 
ie ae j Wi+2CoWi=0. (67) 


The boundary conditions on the surface (r=1) 
are 


i) Wi=0 (68) 
a a? 2rd. z 

ii) E ay me |= 0 (69) 
~ a oe oo 2 
ili) Fas as | Wi=0. (70) 


Without difficulties the solution of (67) is ob- 
tained: 
ie Js/2(air) 
Wie Age ee 
where a, is the minimum value of a;’s which 
satisfy Js/2(ai)=0. 


(71) 


6 
Thus Co becomes C= 3 —4115.48. 


According to CHANDRASEKHAR C> is 3091.4 in 
case of the free surface and 8047.1 in case of 
the rigid one. Therefore, in this respect too, 
it is seen that (63) represents a condition 
between ‘rigid’ and ‘free’, From (66) and (71) 
ru, becomes 


(72) 
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ry = Ahan Vi(0.9) . (73) shown that the above-mentioned direction of 
vr motion is decided in regard to the rotating 
Explicit value of A: in (73) is still unknown @%!S- 
and is determined in the third stage [&]. In From —2rm=p'r1 (74) 
addition to A, there are in fact two more : { 2A Js/o( air) 
arbitrary constants a, bi, contained in Yi(0,g), ™ #8 derived as m= ae art AICS 9). 
tor Yi(0, ¢)=cos 6+a: sin 6 sin y+: sin 6 cos¢. (75) 


But these will vanish when the axis 0=0 is 
chosen artificially to be parallel to the flow >, : A 
at the origin. The situation, however, changes iy general, Ais fapcesse cai ihe a igllowing 


The procedure to obtain v and w is as follows. 


if the sphere is rotating. Later on, ‘it is form; 
ru Was =F EG men) Wrz 0 
FE OW nim r Ln OWmn 
Aca 7 tO “_ sinf a9 (76) 
1 9 Wr OWn.m 
yw See ee eee = Sab LET 
Fam sin 6 0g 7 Lnm 00 


where F, G and L are unknown functions and Wam=r"Ynnl(O, ¢). First of all, L=0, as 
easily seen from (76) and (51). Then the three unknowns are reduced to two (F and G). 
By an appropriate choice of a harmonic Wr,m, 


these two are determined from (73); (F+rG)r=A1 ae (77) 

and from the equation of continuity; - & +r . +4G=0 (78) 

as Ga Aller) a (79) 

The final forms of v1 and wi: are n=Ay sniear) = ms halos) ee 2 (80) 
W.=0. (81) pee 


The general pattern of motion is shown in Fig. 7. 
[€2] The second-order equation is 


2(u,-grad =U (82) 


=a . 
where 21-grad 7: is 


uepeeel d {Leelee} _ — 28" |pcos 0) 
dr 


3a,2 y3/2 Vr 3/2 
21 Az J3/2(air) 1 eared a 83) 
rt 3a,2 y3/2 ES dy idle 3/2 ( 


fy 
The complete solution to the second-order equations 1s - 
Fig. 7. General pattern of motion 


tty = Ag ft cos 8 +0.12A2 fi Dole UIT 


y3l2 


(aur) ar Js/ i ay : Is72(Bir) 
ma—A| Aa — 9 ae in 0—A,?| 0.09 yal 
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A Jsj2(B17) | we (84) 
+0.0372- 7 Seles } fsin 20+ 


elie 
W2=0 
2A. Js/2(@1r) ee 


a? Vr 
0.0043. 4 eee 
. L V y ’ 


Jip(ar) 
Vr 


T= 


s6+0.011A:? PtP) py8)—0.0067 As 


where J5/a(Bn)=0, (2=1, 2:--). 
[€]: The third-order equation is 


V'(rus)—CoL2(rus)+ Ope, +i, -grad t2+u, grad t:]=0. (85) 
J wen Y: (8, 9) 


secular response. Therefore, its coefficient must vanish. This leads to the determination of 
A: which, so far, has been unknown: 


In (85), the inhomogeneous term which has the form of will produce a 


Proceeding as before, we get a solution to the third-order equations: 
us=A A, HEE p, (6)— 0.75 2A) pg) 42.0.4, 2282) po) + 4g AD) Dye) 4... 
2A ) (zr) Jip2(2zr) 
ran Se See P,(0)—0.11A2 ae —0.063As 
—{o.740 sisal) ype ale a \pyo)+0, 18.4,/°” a) P2(6) 
7747) 
+{4 gir) ta +1.6 Ens) +9, pg Hiatt) 4 9y7/tiNrs Poy Pye) +. 
Tsyal ar) ay Ts)2(@ir) OP\(A) —.75 Jsp(@ar) d (20) OP,(@) 
pds | Gan. ek iB 00 aia a yr a0 
Jsi2(81r) “i d Isj2( (fir) OP3(8) 
nh ae ye 3 ar pia a0 
Jsi(r17) Yr a ae) OP3(8) a panty 
+1675] aT + del aan 90 + ‘ w3 =0 (87) 
where 7: are the roots of Jzjo: Jzja(yi)=0, (=1, 2---). 
§8. In case of a rotating sphere, calculations (=-eadae ‘C = Opa on (88) 
become more complicated, as another param- 2 


eter © enters which is indicative of rota- 
tion. The physical quantities are expanded 
in double powers of €andQ. The followings 
are concerned only with the first order of & 
and 2. The chief result obtained is that the 
upwelling site on the surface is decided 
uniquely in relation to the rotating axis of 
the sphere. 

In this case the equation of motion is, 
instead of (51), 


where 0=2oR?/yv, w is the angular velocity 


of rotation and oss (cos 6, —sin @, 0) is the unit 
vector of rotation. 

First of all, we expand C, u, and c in the 
following form, 


C= el + €) 3, Cn" 


me ef St Se (89) 
“= SHS AOM A. , T= Ss By @ Laas ies 
n=0 n=1 M 


m=1 m=0 


A Perturbation Method for Studying Thermal Convection Problems 13 


and then insert (89) into the governing equa- 
tions (49), (50) and (88). The coefficients of 
€Q, become as follows, 


div ws, 1=0 (90) 
—2Fr 1, 1=P'"1, 1 (91) 
a Fa Co 
Pui, 1— grad Wi, i+ 9 ap it et cae of 
=> 
+i, 1X e=0, (92) 


= 
where uw, 0 and 11,0 are respectively previous 


ua and v1 of §7. 
From (90), (91) and (92), 
V°(rt1,1)— CoL*(rur,1) 


(Es => = => 
= “9 L*V*t1,0—PXr-rot rot (u1,0Xe)}=0. (93) 
Hitherto, the direction of the motion iu at the 
centre of the sphere has been undetermined. 
But in case of rotation, it must be decided in 
relation to the axis of rotation. Therefore 


iyo is not given by (73), (80) and (81) as 


before, but 
U1,0 (Fijot+7?Gi,o)7¥i(0-¢)' 
oe 
> F; 0 eae 
us, Mn | = ye? (94) 
oY 
W1, i 
cP Fi "sin 600° 
where Yi=acos0+bsin 6 sin g+c sin 9 cos ¢ 


and a, 6, and c are to be decided. 
Inserting (94) into (93) we get 
p*(rus, 1)— Col? (71,1) 
Jsplaur) { 2C1— 
sal Vr ea 


+c sin @ cos v) +- 


acos 9+ sin @ sin Q 
a? 


on sin 8 cos g 


—csin 6 sin is (95) 


The inhomogeneous terms must vanish iden- 
tically. ipberetore: 


Saree 
Cia=0, oe =0, ake 5 O=0 (96) 
i.e., C:=0, 6=0, c=0. (97) 


This means that the fluid motion U1,0 is paral- 


lel to the axis of rotation at the centre of the 
sphere. 


a 
U1,1 and t1,; are 


Jsj2(@r) 


id eae), al (41,1 CoS 9+b1,1 sin @ sin g 
+¢1,1 Sin 6 cos Y) 


2 Jsj(aur) 
a? y3/2 


Ci,1= (a1,1 COS 8 +b1,1 sin 6 sin g 


+¢1,1 Sin 0 cos ¢) 


Qayy a 5/2, Q\r a F 
ree ha 1 Jenakaar Kain sin 0 


+bi,1 cos @ sin g+¢1,1 Cos 8 cos ¢) 
{[Jsplaur) ar Jsj2(ar)) 


W151 = 4 


| 3/2 2 yil2 pus COS 


—c1,1 sin y)—0.43 Ele) sin20, (98) 


where @,1, 01,1, and c¢i,1, are all unknowns 


which are determined in a later stage. The 
3 5/2, Q1r : . 
general pattern of motion — LO sin 29 is 


axis 
of . 
rotation 


Fig. 8. General pattern of motion 
_ Isjaaar) 


3/2 


sin 29. 


shown in Fig. 8. It is already indicative of 
influence of the rotation. 
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§9. In the remainder of this article, a study 
will be made of the steady thermal convec- 
tion in a two-dimensional fluid layer when it 
is heated uniformly from below under a 
simultaneous constraint of non-uniform tem- 
perature on its upper surface (inhomogeneous 
boundary conditions). The fundamental equa- 
tions are (9) and (10). The boundary condi- 
tions are also (11) except r=0, (y=0). 

As for the non-uniform temperature distribu- 
tion on the upper surface, it is assumed to 
have a periodicity of wave length L, where 
L=nl, | is 2/2, the wave length of spon- 
taneous cells, and m is an integer. Under the 
assumption, it is expected that our model 
(Fig. 9) having a wave length Z corresponds 


d 


le 
ROI RPemIeo.. 
to 


Fig. 9. Geometry of the fluid layer under 
investigation. 


to a spherical shell. In this study m is taken 


as 4. The non-uniformity is expressed as 
follows; 
‘h—\0)5 (Gea) 
res 8V 3 On COS Fonte On), be 
n=l 2a 8Y 2 ws 


(99) 
where 0,’s are arbitray. 
As seen from (99), the non-uniformity is sole- 
ly borne on t; for the sake of mathematical 
technique. 


§10. In this section, calculations are carried 
out as in § 3. 


[ef 


The equations are Or, 


Vidi= 


a 
ay 0 


(100) 


The boundary conditions are 


Ogi Op: Ody 


=o pa a ae to H=0-G=0, 1), 


(101) 


Ox eK’ 


The solution is 


gi =A sin zy sin ——= = ; ==(¢+21) (102) 


== (~+21).(103) 


2 es n COS 7-5 
Bn sin zy ay 


It is to be noted here that besides Ai, an 
unknown phase, %, enters into ¢. and 1. 
This comes from the consideration that the 
site of spontaneous convection cells is to be 
decided according to the surface temperature 
disturbance. In fact, 2: is determined later. 
[€2]: The equations are 


"1S 


OT: 
4 ood eae 
Vid2=Ro x (104) 
Og OT1 Og, Or, Ob. Ps , 
oy Bp On Oy hn Game ae 
The boundary conditions are 
Ops _, Obs Od ee oe 
Ox a Oy? Ox 2 =0, tT2=0, (y=0,1). 
The solution is 
d2= As Sin zy sin [== 7 > = (+a) (106) 
T= ppt F i) Waar wt 
2 . 
-- a sin 2zy. (107) 


In this stage too, a new unknown phase z 
appears. 


&3]: 0 
le | pA Siagee (108) 
Ox Ox 
Od, OT: Ode Or — Od, OT: 
Oy Ox Oy Ox Ox Oy 
Ode Or1 Ow 
[ a Oe Oy Fd Aa (109) 


The boundary conditions are from (99) and 
(11) 


Ons Obs Obs 
Ox’ Oy? Ox? =o, w= 0D), 
Ta Os) (110) 
8Y 2° 
TA og ey On COS wr n(z+On), (y=1). 
(111) 


In the first place, (108) and (109) are cross- 
differentiated and then subtracted in order to 
get an equation governing qs. 
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hi Ors 
Ea Er; 
pos 4 ae 
7c? A? F : 
5 AR| of =1 fsin my sin ¥e 5 -(a+211) 
—1352A;,' sin 3zy sin — rae 5 —= (a+a1). (112) 


Considering (111) we seek a solution of the 
form 


ps= Sal) sing n+ On) 


+ 0.000164.A,3 sin zy sin —— 


=(%-+a1), (113) 


a es 
where the last term in (113) is the particular 
integral of —1352.A,3 sin zy sin eee (a+2a1) 


of (112). 


As the following manipulation is somewhat 
complicated, it will help us to rewrite the 
boundary conditions using fn’s: 


ed =U “ty=0, 1) (114) 
r= 0) (y=0) (115) 
FD (y=1) (116) 


Now the problem is concerned with the func- 
tions fn’s. Of all these fn’s, fs is the most 
difficult one to deal with. 

i) fs. From (112) and (113) it is seen that 
(117) 


m=O, 
e., the phase, 21, has been determined. 
The equation governing /; is 


@ 
la | fan + “Fay 
Sars —1 fsin xy (118) 


The general solution /,(y) of (118) is the sum 


6 
of the homogeneous solution a ahiy) plus a 


particular integral /,(y); 
= 3.70 AiRo | Ay? 
Fi=— aly 


2% 4 COs zY, (119) 
10°73 | 24 1 jy ie 


Faly) = E ashy) 


is 3.70 AiR, A? 
10273 24 
where ai’s have not been determined as yet, 


1.ée., 


— 1| y COS ty, (120) 


and /i(y)’s are the solution to the homogeneous 
equation of (118); 


a 
| aes [m+ “Yh y)=0, (121) 
(u,-+++he)=(sin zy, cos zy, e94+1,62 Oy 
DAE ABR CONFESS EET Cam Let Se (122) 


Next we have to make (120) satisfy the 
boundary conditions (114), (115) and (116): 


6 
> ahi(0)=0, 
i=2 


—3.70AiRo Eo 
len? 124 


S 
a 
ie 
l 
= 


(123) 


iMe> iMe 


=3 204 Ry As i 
2 li B= Bab 


aihi(0)=0, 


6 Sister 
> adel) = 
—4R ods , 


where the summation on the left-hand side is, 
as noted above, from 2to6. This is because 
hi=sin zy vanishes on the boundaries and 
need not be written down. In (123) there are 
six linear equations while the number of un- 
knowns ai:’s are only five. Therefore if (123) 
is to be consistent, all of the equations can- 
not be linearly independent. There must be 
six of such &’s as, when multiplied by the 
equations (123) respectively in order, and 
summed up, will make the coefficients of ai’s 
vanish: 


E,hi(0) + Eohi(1) + Eshu(0) + E,hi(1) + £5h(0) 
+E hi (1)=0, 6). (124) 


From (124), an are determined 


1 


—3.70Ai Rox E | 
24 


(W=ee . 


10 4 
Gitte a a $4 — Tana 68 Tar 
4 
f= Fa (125) 


where &; is taken as 1. 

One more thing must be added. If (123) is 
consistent at all, these &:’s must also satisfy 
the tollowing relation: 
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3.70 Ai Ro ae aye 

—4R6.€6= (126) 
which is derived from the inner product of 
the right-hand side of (123) and &’s. (126) 


contains A: so far unknown. This means 
that (126) is the equation determining Aj. 
Inserting (125) into (126), the final relation 
determining A; is obtained: 


of 4 


27 
54 (127) 


The relation (127) is represented in Fig. 10. 
In some cases, for a given 0, there may be 
three A:i’s which formally satisfy (127), e.g., 
and A,’’ 


Ai, Ai’ in Fig. 10. Of these three 


Fig. 10. Curve for deciding the amplitude of a 
spontaneous cell. 


A,’s, only A: is the one we want, for with 
continuous increase of 04, Ai’ and Ay’ sud- 
denly cease to exist. This means that Aj’ 
and Ai” have no physical counterparts. 

The procedure to get to the relation (127) 
was somewhat circuitous. In another simple 
way, (127) is also obtainable. As in case of 
the Grreen’s integral theorem, we first con- 
sider the homogeneous equation of (118) ; 


OKAI 
a fat Nee 
& [sin nil dy? —= }sin cy 
= (5) sin ry=0 (128) 
with a solution, sinzy put into it. Then 


multiplying (128) by fs(y) and (118) by sin zy 
respectively, integrating them both from 0 to 
1, and finally subtracting one from the other, 
we get 


Rez G [Siem] Sin ae 


7? Ro al -1| 


4 
This is the same as (127), as easily seen on 
integrating the left-hand side of (129) by parts. 
Although the present method is easier, we 
cannot fully rely on it, for from (129) the 
functional form of /s(y) cannot be obtained. 
il) fr (n*4). The governing equation is 


A? 


3A (129) 


P°faly) sin n(2+On) 


ve 


27 
— Fat fal) S> os 575 M(e+On)=0. (130) 


1a 


The boundary conditions are 


fay =0, faly=0, (y=0, 1), fay) =0, (y=0), 


(131) 
frly)=—nRodn, (y=1). (132) 
We seek a solution in the form 
3 
Sry) = x ctet ry , (133) 


where ci*’s are unknown constants and /pi’s 
are the roots of 


5 a 2in 2 
(p 39 ur 198 ” >=0 (134) 
her. 

rf nen 323 (n? 1/3 _ nx? 
pint a (a) Soe Ba 


3772 2\ 1/3 = ; 
+ (=) I+Y 31)=(@+ri)?, ps*=(q—ri)?. 
(135) 


Taking the boundary conditions (131) and (132) 
into account, we get for fn 
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pris NRoOn sin sy j Tt 
n “ark(s*-hq?-+ 7°)? —4g?y?) {(8+¢ +r) 4g) eee P — 0.000986 A.’ sin 3zy cos Mets 
BL 08 se a a(y+1) 
V 3 (e+ e-"0—2 cos 2r e cos( ry—1+ = +0.00623 A: sin zy cos —- 7a —=(x%+0,) 
es T << i 
+ e-a9 cos (y= ee 5) =—gaiy ot) cos( ry igh: 5) rie 2 Ongn(y) COS ——— 8 aT n(x+On), (138) 


—e-9"-Deos{ ry-F1— 5)| ] (n24). (136) 


To sum up, ¢; becomes as follows; 


$3= A3 sin zy sin S== (#+2:3) 


a 


+0.000164A,8 sin 3zy sin—== (#+ 64) 


oe 


On n > n 

+ >) Fnly) sing 5 5ue+O ys (137) 
where fn’s are shown in Fig. 11. Of these 
Fn’s, fs and f; are larger in value than other 


Pig. 12. 77s) in (60): 


where gn’s are shown in Fig. 12. Among 
these, g3 and gs are large by far. 


[€4]: The equations are 


Ot.s Ot: 


dhe 
Vi du= are Di Dre 5 


Od, OTs Ode Ort2 Od: Ot1 = Od: Ot: 
Oy Ox Oy Ox Oy Ox Ox Oy 
Obs OT. Obs Or1 Ode a, 
Lah Soe Oy OPA Die 
The boundary conditions are homogeneous. 
From (139) and (140) we get an equation 
governing ¢, and then the relation determin- 


(139) 


ing Ao: 
Fig. 11. f;’s in (48). Asx? [42-1 -o. yr eth (141) 
Ffn’s. This is due to the fact that, in case of /s, : : 
and fs, sins in (136) is nearly zero. As to the convergency, it is seen that the 
ts is obtained from the auxiliary equation terms of [&*] are negligible so long as A,~8 
(108) or (109). and &*<1. But the range of € is narrowed ' 
A 1A still more in the stage of [&°], as inhomogene- 
t3=— co 7 9 (G4— a4) sin 2rry ous terms containing 63; and 6; produce a 
particular integral of a large value. There- 
V 2 As fore in dealing with the 03 and 6d, disturb- 
3% se 7 COs V Ty ete) ances, enough consideration must be paid. 
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§11. It is of interest to see how the patterns 
of convection are changed by the inhomogeneity 
constrained on the surface. This will be seen 
from the following several examples. 

i) The temperature difference between the 
upper and lower surfaces is 10% in excess of 
the critical value. The surface temperature 


disturbance is in the form cos a x and its 


amplitude is 10% of the temperature difference 
(6,=1.75). From (127), A=—7.35. Isother- 
mals and stream lines in this case are shown 
in Fig. 13. The site of spontaneous cells is 


Sabena s Sse 


2020 2e=42 
Fig. 13. Isothermals and stream lines when the 
temperature difference between the upper and 
lower surfaces is 10% in excess of the critical 
value. The surface temperature disturbance is 


v1 
of a form cos /2 
the temperature difference. 


« and its amplitude is 10% of 


Temp. diff 


OKAI 


decided entirely according to the inhomogeneity 
on the surface. The upwelling site is where 
the upper surface temperature is the highest 
and the sinking site is where it is the lowest. 
But, the pattern itself is almost the same as 
that with no inhomogeneity on the surface. 
It is to be noted here that the inhomogeneous 
temperature disturbance having the critical 


wave length, 2)/ 2 in this case, decides once 
and for all the site of spontaneous convection 
cells, however small the disturbance may be. 

ii) The temperature difference between the 
upper and the lower surfaces is 10% in excess 
of the critical value. The surface tempera- 
ture disturbance is 


1 al 713 
0.17 az; [cos eget 5 oS e| E 


Isothermals and stream lines in this case are 
shown in Fig. 14. It is seen from Fig. 14 
that surface disturbances having much larger 
wave length than the critical one play very 
little part. Neither high nor low temperature 
sites on the surface seem to have considerable 
influence on a convective pattern. Much less 
can it be said from Fig. 14 that convective 
currents are upward where the temperature 
is low on the upper surface and downward 
where it is high. From this fact upward con- 


Cold 


— SG SS = 
SRO) LESS LEE 
: SEZ i SEZ | 


Fig. 14. Isothermals and stream lines when the temperature difference between the upper 


and lower,surfaces is 10%, in excess of the critical value. The surface temperature 


v1 


. , 1 T ; 
disturbance is 0.1 Pain| cos P52 3008 75° | in other words, the amplitude for the 


wave length 8\/ 2 is 10% of the temperature difference, while that for the wave length 
2/ 2 is half as large. 


A Perturbation Method for Studying Thermal Convection Problems 19 


' 48% 


Fig. 15. Stream lines when the temperature difference between the upper and lower sur- 


faces is only 0.3% in excess of the critical value. 
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vective currents may not always be expected 
under the oceans where the temperature is 
comparatively low compared with under the 
continents. 

iii) Surface disturbances having wave 
_lengths close to the critical one have much 
influence in determining the general feature 
of fluid motion. In §10, these are terms con- 
cerned with 63 and 0d;. 

Shown in Fig. 15 are stream lines when the 
temperature difference between the lower and 
upper surfaces is only 0.3% in excess of the 
critical value and the surface temperature 
disturbance is 

37 1 ae 

A 2 z+ o (1) cos get 

In this case, as the surface disturbance is 
very small, its influence is not seen very 
clearly. But if we compare Fig. 15 with Fig. 
14, it is seen that the disturbance of 03 is 
very influential. 

Fig. 16 shows the root mean _ square 
amplitude (in km) of a spherical harmonic of 
order 7 in the topographic analysis of the 
earth’s surface after Vening Metnesz (1951). 
CHANDRASEKHAR (1953) showed that the convec- 
tive patterns of harmonics of orders three 
and four set in first when the depth of spheri- 
cal shell is about half the radius of the sphere. 
Then he indicated the bearing of this result 
on the problem of convection in the earth’s 
mantle and of the interpretation of the earth’s 
topographic features. In this present article, 
it was shown that surface disturbances having 


~ 0.003 Taisy COs 


The surface temperature disturbance 


x, in other words, the amplitude for the wave length 


is 0.3% of the temperature difference, while that for the wave length 2,/2 is in- 


much larger or smaller wave lengths than 
the critical one play very little part in deter- 
mining the general feature of fluid motion, 
while those having wave lengths next to the 
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Fig. 16. Root mean square amplitude (in km) of 
a spherical harmonic of order J in the topo- 
graphic analysis of the earth’s surface after 
VENING MEINESZ. 


critical one are quite effective in this. This 
conclusion is of some interest as it is seen 
from Fig. 16 that the amplitude of a spheri- 
cal harmonic of order 5 is considerably large 
compared with that of order 6. 


In conclusion, the writer wishes to express 
his sincere thanks to Prof. C. Tsusor and Dr. 
H. Taxeucut of Tokyo University for their 
encouragement and support throughout the 
course of this study. His thanks are also due 
to Prof. S. CHANDRASEKHAR Of Chicago Univer- 
sity for calling his attention to Matxus and 
VERONIS’ study. 
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Abstract 


At first, some examples are calculated of the dispersion curve, the amplitude func-— : 
tion and the amplitude distribution for the first order of LovE-waves. These calcula- 
— tions remained for cases (A) and (B) in eq. (57) in the previous paper®). 

j Next, for cases (C) and (D) in eq. (57), the other examples for the zeroth and the 
first orders of LovE-waves are presented. : 

Lastly, some physical considerations as to these calculations are given. 


13. Numerical calculations of the dispersion curve for 
L, in cases (A) and (B) ; 
x The dispersion curve for the zeroth order of ‘LoveE-waves, noted by bas 
were already calculated for cases (A) and (B) of (57) in Sec. 9. The rigidity 5 
_ ratios adopted in the numerical calculation were (a) p2/i:=4 and ps/pu2=30, 
(b) appa—30 and Lin} fia —= 30, (c) pal 1/16 and 183 / foe = 30 and (d ) jia/ Pa=et/o 
and j3/j2=30. 
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The dispersion curve for the first order of Love-waves, noted by L,, has 
been calculated in this section, adopting the same rigidity ratios as before, 
The results are shown in Figs. 7(a) to 7(d) for various values of the parameter 
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c=U=Vs 


Tv 7H 
Fig. 7 (dj. Ly, 
Fig. 7. Dispersion curves for L, in cases (A) and (B). 
(a) Bo/Py = 4 and p,/p, = 30, (Db) Hy/Hy = 30 and ps/u, = 30 
(C) Bo/fy = 1/16 and ps/p, = 30, (d) fy/my = 1/3 and p,/H, = 30 


H,/H,. 

_ Whereas L, has no cutoff period, L, has a finite cutoff period at which 
phase- and group-velocities must coincide with the velocity of S-wave in 
the lowest layer. Periods of L, must be confined between the cutoff and zero. 


a 


O taint ats ees tan hal Weel L Nl eae oe 
OAl 1.0 me) ey aye 
Fig. 8. The period corresponding to the second minimum group velocity 
for L, in cases (A) and (B). T’ = T/(H,/v,+ H,/2,). 
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Group velocity has two minima as in the case of the zeroth order; one is 
near T/(H,/v,) =4/3 and another is near T /(H ,/0,+H2/v2)=4/3, for the 
rigidity ratios (a) and (b). The latter relation will be clear from Fig. 8 where 
[lg/ J42 1S kept at 30. 
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Fig. 9. Amplitude functions for L, in cases (A) and (B). 
T’ = T/(Hy/0,+ Hef 09). 


14. Numerical calculations of the amplitude fuuction 
for L, in cases (A) and (B) 

The amplitude function 27A(&) in (45) has been obtained by using the 
dispersion curve illustrated in Fig. 7. The results are shown in Figs. 9 (a) 
to 9 (d) as functions of period Tv,/H,. 

These amplitude functions may have two maxima in each case. It must 
be noticed that the first maximum in cases (c) and (d) occurs respectively 
near the period of the maximum group velocity, whereas the other maxima, 
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Fig. 10. The period corresponding to the second maximum of 
the amplitude function for L, in cases (A) and (B). 
T’=T /(Hy/0,4+ H2/v2)- 
(a) : (4, 30), (b) : (30, 30), (c) : (1/16, 30), (d) : (1/3, 30) 
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the first maximum in cases (a) and (b) and the second maximum in all cases, 
occur respectively near the period of the minimum group velocity. 

Fig. 10 shows the relation between the period of the second maximum of 
the amplitude function and H,/H,, being parameters (ps/ji, Piaf pin) np Lt. 16 
clear that each curve will approach to the period indicated by T/(H,/21 
+H,/v,)=4/3 with the increase of H,/H,. 

Now it has been found that either L, or L,; must obey the wave-length law 
in general, 


T/(H,/01) or T/(E,/t+Ha/o) =4/Qn+1), ("for be) (4) 


4 
at least qualitatively. The thickness of the layers may have larger influence on 
L, than on L,, because the former has shorter wave-length than the latter. 
Comparing Fig. 10 with Fig. 5, the asymptotic value of (64) seems to be easily 
attained by L, rather than by Ly. But the ratios of deviation from the 
asymptotic value are nearly equal to each other for the same value of H,/H,. 

This situation may also be seen in the relation between the sharpness 
of the maximum of the amplitude function and the period. The higher the 
order of the waves, the narrower the band of the period having considerably 
large amplitude. 

It will be physically expected that the amplitude function may arrive at 
the maximum near the period of the mimimum group velocity. Analytically, 
this phenomenon may be easily recognized by the factor (v,/U—v,/c) in (45). 
On the other hand, it will seem strange that the amplitude function may arrive 
at the first maximum near the period of the maximum group velocity in cases 
(c) and (d). Numerically, the present authors have found that the factor 
{1—(v,/c)?}-1 in (45) plays the greatest role in these cases. 

The latter phenomenon must be accompanied by the “low velocity 
layer”. Considerable energies will escape from the most superficial layer 
to the low velocity one. Thus the wave of short length cannot keep the law of 
(64), but ignores that law. Displacement of the wave of short length has 
the factor cosh #,z as shown in (iii) of page 126. If c¢ coincides with v, 
cosh #,z becomes unity for each depth, resulting in large displacement on the 


surface of the earth for a certain amplitude on the lower boundary of the 
most superficial layer. 


15. Numerical calculations of the amplitude distribution 
for L, in cases (A) and (B) 


Figs. 11(a) to 11(d) illustrate the results calculated for several periods 
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Fig. 11. Amplitude distributions for L, in cases (A) and (B). 
T’ = T/(A,/v,+ He/v). 


: , ) and 11( (b), in the most superficial eee for” ‘ ae 
Oxenk the other hand - in Figs. I1(c) and 11(d), the nodes migrate 
ards A the low Os layer and go into the upper layer with mich nes > 


: > Ketaal amplitude ceenpatinn does not always show jancean configura- 
tion forming 3/4 wave-lengths, whereas (64) is satisfied fairly well at the 
maximum of the amplitude function. 
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a 16. Numerical calculations of the dispersion curve 
for L, and L, in cases (C) and (D) 
pr Case (C) from (57) may take the type (iii) in (26). The dispersion curve 
. has been calculated by the characteristic equation (iti) in (27). 
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Fig. 12 illustrates the first two modes, L, and L,, in the two cases (a) 
ps/ pa =1/80, bs/pe=4 and (b) pafj4=1/30, jee] y= 10; the parameter being 
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Fig. 12. Dispersion curves for L, and L, in cases (C). 
(a) g/t = 1/30 aud jug/p, = 4, (b)  2/, = 1/30 and jy/p14 = 10, 
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Fig. 13. Dispersion curves for Ly and L, in case (D). 
By/ ly = 30 and fs/Hy = 1/4, (b)  Hy/Hy = 30 and ps/P2 = 1/10. 
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_ Fig. 14, Amplitude functions for Ly and L, in case (C). 


Blac amplitude function is slight when H, /H, is small. However the maximum 
value increases suddenly with increase of H,/H,. 
In case (D) the amplitude function for Lo, as shown in Fig. 15, has the 
- maximum. near Tv,/H,=4 and that for L; near To Hi =4/3. 

_ It may be curious that the amplitude function for L, in case (D) is not 


- influenced by the ratio H,/Hy. 
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18. Numerical calculations of the amplitude distribution wv 
for L, and L, in cases (C) and (D) 


The results are exhibited in Figs. 16 and 17 for several periods indicated — 
by T/(H,/v,) or T/(H,/v,) respectively. 
The amplitude distribution in case (C) may form a type of half wave | 
-length within the second layer, whereas that in case (D) may form a type of we, ~ 
quarter wave-length in the first layer. ; a 
It must be noted in case (D) that the amplitude in the lowest layer. be ; “aoe 
comes remarkably small, if j3/j1. decreases or if H,/H, increases. ; : 
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19. Conclusions 

It is the principal aim of this paper to investigate the basic property of 
LovE-waves. In order to make the idea concrete, rather large rigidity ratios 
are assumed in the numerical calculation, regardless of actual conditions of 
the earth. These large rigidity ratios will be found only near the surface of 
the earth in seismic prospecting. 

“Quarter wave-length law’ has been often experienced in actual observa- 
tions”’®), It has been expected physically that the amplitude of displacement 
would be large at the phase corresponding to the minimum group velocity. 
However this expection was too abstract to be connected with that law 
actually experienced. 

If a pulse is generated from a source, the observed displacement must 
have the factor {|dU/dw|}—1’2, as is well known theoretically?) Indeed the 
displacement of the stationary phase of group velocity will be large, owing to 
this factor. But the contribution of this facter is not so effective in making 
the amplitude of the stationary phase larger than that of the other phase. 

The present authors have noticed the amplitude function 27A(€) itself 
at periodic motion. They have found that the law of (64) was satisfied by the 
phase corresponding to the maximum of the amplitude function. Due to 
the factor (v;/U—v,/c) contained in 27A(E), it will be recognized moreover 
that the minimum group velocity may give the maximum amplitude. 

(64) recalls to our minds the vibration of sound in a tube, one end closed. 
But it was obscure, with this analogy alone, why the period corresponding 
to the minimum group velocity would satisfy (64). To tell the truth, a 
quarter wave-length in z-direction is shown by the amplitude distribution 
more completely at Tv,/H,-+0 than at Tv,/H,+4. It must be remarked here 
that the true wave-length in z-direction is 27/7, and Tv, means nothing but 
an apparent wave-length. 

In spite of this physical uncertainty, the acuthors think, the calculated 
result (64) must be interesting. The observed experience and the abstract 
consideration, mentioned at the beginning of this section, have been connected 
well by (64) to each other at least numerically. 

In some cases the maximum of the amplitude function is got also near 
the phase corresponding to the maximum group velocity. This phenomenon 
is physically less clear than that described already. The low velocity 


layer must play a grand role in it. 
Shielding effect ‘Quarter wave-length law” is apt to be satisfied by the large 
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velocity contrast between layers. The nearer the surface ‘of the earth the 
boundary exists, the larger becomes its effect. If p2/j1=4s/p2, for instance, 
the third layer cannot play an important role but the second layer plays the 
greatest part. This phenomenon may be recognized by comparison of Fig. 6 
(a) with Fig. 6 (b). In Fig. 6 (b) very small energies penetrate into even the 
second layer. In this circumstance the boundary between the first and the 
second layers acts as if it were a shielding screen for the lower layers. 

In Fig. 6 (a) the subsurface has two screens; one is the boundary between 
the first and the second layers and the other is that between the second and 
the third layers. The latter screen is more effective than the former in Fig. 
6 (a). Of course, the thicker the layer acting as a screen, the stronger the 
shielding effect. 

Owing to the shielding effect, it will be easy to deduce, by the observation 
of surface waves, the geological construction near the surface of the earth. On 
the other hand, deep constructions cannot be deduced from the observation 
of surface waves, because deep layers must be masked by shallow layers. 

Surface waves generated near the surface must be mainly affected by very 
small constructions near the surface and they will take nearly whole parts of 
“noise” on the ground. 

If the source of the wave lies deep from the surface, only the waves 
having large wave-length must be taken into account, because the waves having 
short wave-length have slight amplitude on the ground in this case. This may 
be understood by Fig. 6, recalling “reciprocal theorem’’!). LoveE-waves having 
as large wave-length as the thickness of the earth’s crust have been observed 
in natural earthquakes*), Natural earthquakes have so deep sources that 
the surface waves having short wave-length which might be affected by small 
construction near the surface cannot be observed on the ground, no matter 
how large a velocity contrast exists near the surface of the earth. 

Now we will conclude that the first boundary which has a large velocity 
contrast under the source plays the most important role in growth of Love 
-waves. To say it in other words, that first boundary alone will be found by 
the observation of LOVE-waves. 
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Transition from Dispersive RAYLEIGH Waves to 
Sound Waves in a Layer over a Half 


Space Absolutely Rigid 
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Abstract 


Dispersive RAYLEIGH waves in a layer over a half space absolutely rigid were 
investigated pretty thoroughly by GrEsel). 

He ignored M')-type of the waves, assuming that this type had much smaller 
amplitude than that of M“)-type. However the present paper suggests that his 
assumption had no theoretical foundation. 

Moreover he did not refer to any higher order of M-waves. Indeed he treated 
media having Polsson’s ratios larger than 0.25, but the character of M-waves was 
obscure at the limit when Potsson’s ratio arrived at 0.50. 

It is the main purpose of the present paper to clear up the last defect just men- 
tioned. 

It must be noted that notations and expressions in this paper will follow those 
employed and obtained in the previous paper by the present author?). 


1. Superficial waves in a liquid layer 


Taking space-coordinates as shown in Fig. 1 and assuming a line source 
at z=E, displacement-potential of superficial waves can be written as 

( be tee : 
Ee ele . ) Sina E sinaz, (1.1) 


[pb |\Y-mMo= — 2A { U 


/ 


air 
me) mx > X 
E 
Me 
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aes aap 
absolutely 
rigid 
v 
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Fig. 1. Space-coordinates under consideration. 
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being the same expression as that obtained in the case of a plate®?. 
In (1.1) the characteristic equation is given by 


M“) —0 thatis cosa H =0 (1.2) 
or eH = (OL 9/2 = 01d eee (1.3) 


2. Dispersive RAYLEIGH waves in a solid layer 


In a solid layer, as well known, M@) and M®)-waves must exist. The 
displacement-potential of them can be expressed, by means of a method similar 
to that used in a solid plate, as follows. 

(i) c2u,>0,; =a, B= B. 
ug 1 1 2 Poe > 
[p]M-o = 2H ( ai ; ) [cos a (¢—E) {(1—AA’) sin aH cos BH 
+(1+AA’)cos aH sin BH}+sina(z+E) {—(A—A’)cosaH cos BH 


+(A+A’) sin aH sin BH}+cos a (z+E) {(A—A’) sin aH cos BH 
+(A4+A’) cos aH sin BH}]/N, 

[Wln-o = 4 5 oy e +.) [B {(14+.4’) cos a (E—H) cos B (z—H) (2.1) 
—(1—A’) sin a (E—H) sin B (ez -H)}+B’ {(1+A) cos aE cos Bz 


—(1—A) sin aE sin Bz}]/N , 
N = {(1—AA’) + (B/a) (1+.4A’)} sin aH cos BH +{(1+AA’) 
+(B/a) (1—AA’)} cos aH sin BH , 
in which the characteristic equation is given by 
M=0 that is sin?}(a+8)H—AA’ sin* } (a—f) H=}(1—AA’—BC’). 
(2.2) 
A’ in (2.1) and (2.2) means PP-reflecting coefficient as to displacement 
-potential on the lower boundary in a superficial layer. 
Letting rigidity of the layer tend to zero, A must approach to -1 and B 
to zero, as already mentioned in the case of a plate. On the other hand, A’ 
will approach to 1 at the same time, whichever value phase-velocity may have, 
as illustrated in Fig. 2. B’ as well as C’ will approach to zero again, being 


B"?—(a/B) (I-A). In this case, [yr]m-o in (2.1) must be zero and (2.2) will 
be reduced to 


{(tan oH /2—1) (tan BH/2+1)} {(tanaH/2+1) (tan BH/2—1)}=0, (2.3) 


that is 


aH = (21+1)27/2 or BH =(2m41) 7/2, (2.4) 


= O/ Ve 


Fig. 2. PP-reflecting coefficient on the lower boundary in a layer. 


Be where J and m represent positive integers, indicating orders of these curves. 

i The former of (2.4) coincide with (1.3) and one sees that [fl] 9 in (2.1) will — 
fic coincide with that in (1.1). 

Jaa . On the other hand, when the latter of (2.4) is satisfied, [JM in (2.1) must 
es ‘be zero. 

¥ fe, Thus eae rigidity of the layer tend to zero, branches corresponding to 
ea (tan «H/2-1) (tan aH/2+1)=0 alone can survive, otherwise amplitudes of — 
displacement-potentials become zero. The surviving two branches coincide 


Bs 

et respectively with the branches of the even and odd orders of superficial waves 
eo in a liquid layer. 

i) -v,2e2B,; a=—ia, p=s 

AS iy =P 


[plu-o= 7% (4. — 1.) [cosh « (eB) {et sin (BH—-8-8) 


fe-4¥ sin (BH +8+8 )}—{et#+-™) sin (BH 8—8 ) 
fe*@+E-) sin (BH +8—8')}]/N, 


wate (4-1 Y 28 rinz5pn[a0em snip tA) Las 


as ATAU eC 

: —e-E-H) sin (B (2 —H)+8’)} +(sin 2 5’)1? {e#” sin (Bz— 6) 
vie —e-*F sin (Bz+8)}]/N, 

es 


NN Se sin (B H-8—8') + e-*” sin (B H+5+8’) 
é +(B/a) {e cos (B H—S—8')—e-* cos (BH +5+8')} , 
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in which the characteristic equation is given by 
M = 0 that is 
{es#/2 cos} (BH—0—0') 2-4 cosas (BH+6+8’)} 
- {et#/2 sin (8B H—8—8 ) +e-*4/? sin $ (BH+6+8')} 
= sin (5+6’) + {sin? (6 +6’)—sin? (6—8’) }*/?. (2.6) 
6’ in (2.5) and (2.6) means phase-lag of displacement-potential at PP 
-reflection on the lower boundary of the layer. Letting rigidity of the layer 
tend to zero, 6 must approach to zero, as already mentioned in the case of a 
plate. On the other hand, 6’ will approach to 72/2 for all values of phase 


-velocity between v, and v, as illustrated in Fig. 3. In this case, [W]m-» in 
(2.5) must be zero and (2.6) will be reduced to 


cosh 4H cosBH =0. (2.7) 


Because cosh 4H cannot be zero, cos 8H must always be zero. Owing to 


bE d= 
OF oe 
es 
~ i ya 
9 J 7 
es f 
ye rae 
| / é , 
/ 
/ 


0.0 Ol o2 O03 O4 OS Os O7 08 Og Ke) 


+ p 


Fig. 3. Phase-lag at PP-reflection on the lower boundary. 


—E){(1—AA ) sinh @H cosh @H 


+(1+A4’) cosh aH sinh 8H} — sinh a (2+ £) 

_ + {(A—A’) cosh &H cosh BH +(A +A’) sinh 4H sinh 6H} 
+ cosh 4 (2+£) {(A —A’) sinh af cosh 6 H 
_-+(A+A ) cosh 4H sinh BH}]/N, 


2 nae er o> fot 1 z (2.8) 
Whee 6 SE (4-1) Bta+4 ) cosh « (E—H) cosh p (2H) 3 
ae. +(1—A’) sinh ¢ (E—H) sinh 8 (z—H)}+B’ {(1+A)coshaE | 
Bp - cosh z+(1—A) sinh &£ sinh AYN, Bh 
ae. N =—i[{(1—AA’)+(6/8)(1+AA’)} sinh 6H cosh BH toe 


— [(14.4A )+(8/&) (1-AA’) cosh &H sinh BH , 


baer which the characteristic equation is given by 


C/Vp 
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Fig. 4. PP-reflecting coefficient on the lower boundary. 
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M= 0 that is 
sinh? } (4+) H—AA’ sinh? 3 (a—f) H =4(BC’+AA’—1) =P. (2.9) 
A’ in (2.8) and (2.9) means again PP-reflecting coefficient on the lower 
boundary of the layer. 
Letting rigidity of the layer tend to zero, A must be confined to —1, as 
already mentioned in the case of a plate. On the other hand, 4A’ is to be 


confined also to —1, as illustrated in Fig. 4. Because B and B’ are again zero 
in this case, [yy]m-) in (2.8) must be zero and (2.9) will be reduced to 


cosh (4/2) sinh (BH /2) sinh (4H/2) cosh (BH/2) = 0. (2.10) 


As cosh (aH/2) sinh («H/2) cosh (8H/2) cannot be zero, sinh (BH /2) must 
be zero, with the result that [],;-) in (2.8) is also zero. 

It has been recognized in this section that dispersive RAYLEIGH waves in 
a solid layer should be reduced to superficial waves in a liquid layer at the 
limit when rigidity is taken as zero. In the following section, processes of 
transition will be investigated. 


3. Numerical calculations of the dispersion-curves 
for various POISSON’s ratios 

At the limit when Porsson’s ratio arrives at 0.50, A’ in Fig. 2 must 
coincide with 1 and A with —1 for any value of c/vs. 

Turning attention to A’ and A for general Porsson’s ratio, one sees that 
A’ will always approach again to 1 and A to —1 if c/v,>1. 

Considering the above two circumstances with Fig. 2, it may be expected 
that the nearer o approaches to 0.50, the wider becomes the range satisfying 
(2.3). 

If c/v,=1, one sees A’=A=~—1 for any value of o, with the result that 
the right hand side of (2.2) becomes zero, because BC’=0 in this case. 
Namely one has 


sin «H sin BH = 0, 
being the same as the relation in the case of a plate. 
If c/v,=1, one sees again A’=A =~] for any value of ¢. Thus one has 
sinh &H sin BH = 0. 


Because sinh 4H cannot be zero in this case, dispersion-curves must always 
approach to sinBH=O0 near c=v,. 


As one has seen by now, for any value of as 
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cosaH =0 and cospH =0 (3.1) 
might be the bases with respect to dispersion-curves of RAYLEIGH waves in a 
layer resting on a half space absolutely rigid, if c/v,>1. On the other hand, 
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Chain-lines correspond to aH =(2/+ 1)" /2, while 
dotted lines to BH=(2m+1) 7/2. o=0.48. 
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Fig. 5. 
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sin BH =0 (3.2) 
might be the base near c=v,. 

In Fig. 5 are shown the curves given by (3.1), that is (2.4). Classification 
of M@) and M®) in Fig. 5 follows that in the case of a plate. Odd / and even 
m belong to M®), on the contrary even / and odd m toM®). The curves given 
by cos BH=0 must be shifted to the rightward with increase in o, though the 
curves given by cos aH=0 need not. In Fig. 5, for an example, the curves 
for «=0.48 are exhibited. 

(i) Now /#2=4 (1-AA’-BC’) in (2.2) has been calculated at first and is illus- 
trated in Fig. 6. If c/v,>1, one sees that J is equal to unity for any value of 
o and (2.2) may be resolved into the next two branches : 


0.0. | 
1.0 1.5 2.0 25 


Tie, Wp 


Fig. 6. The relation between / in (2.2) and C/V». 


es sls + ee +(—AA )*/? sin ( : a +) EH 0 for M{), 
cos ( -F Ze E | ou —(—AA)é* sin ( — +.) # = for M(&) (< 3) 


where classification of M“) and M®) depends on the limit when (3.3) will 
coincide with (2.3) in which A=-land A’=1. It will be very easy, asin the 
case of a plate, to obtain dispersion-curves from (3.3). 

If c/v, approaches to 1, / differs from unity and (2.2) cannot be resolved 
into factors. The process of graphical solution for this case, or 


in( Bg) Has oanrene( fg)" oy 


An example of graphical solutions for (2.2). 
=0.48 and ¢/vp=1.5. ; 


is shown by Fig. 7. Owing to intersection of cos sH=0 with cos BH=0, 
systematic classification of M® and M®) is somewhat disturbed, if c/v, is 


Ma not much larger than unity. This disturbance appears in Fig. 7 between : i 
Be.  M,® and M,®). 

% a (ii) (2.6) may be rewritten by 

a 5 sin (BH —5—8’) = e—*44 sin (BH +5+6’) + 2 e~*# (sin 2 6 sin 2 6’)!/2 

ce ; ioe for c/u,=V2. (8.5) 

| ‘ Re An example of graphical solution for (3.5) is exhibited in Fig. 8 where 

ee M,,) and M,,?) appear in turn regularly, because cos aH =0 does not exist | 

; in the present region and no disturbance can occur. 


Mo Mo 
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Sed 4y 


Fig. 8. An example of graphical solution for (3.5). 
o=0.48 and c/vp=0.50. 
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© Fig. 7 alone. In order to avoid this difficulty, 
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(iii) -« (2.9) may be rewritten by rh 
‘ 7 ay f . & B EH iA ese = & ti Bp Eo 1/ 
stay . sinh ( E ‘ 22 ) meas \p+AA sinh ( E ) 5 


from which M, alone will be obtained. ao 
If c becomes smaller than v,, A and A’ will decrease to Become remarkably 

smaller than —1 at once and /? in (2.9) will coincide with AA’. In this case 

(3.6) will be reduced to ; wi coe 


sinh ( : - : cH =1cosh ( . p ) cH d 3.7) aes: 


2 Ly tee zZ 
As / in (2.9) is illustrated in Fig. 9, graphical solution of (3.7) is very easy. . 
Thus the dispersion-curves for «=0.48 have been obtained, as shown by 
full lines in Fig. 10 where the other lines correspond respectively to those in 
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Fig. 10. Dispersion-curves for o=0.48. 
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Fig. 11. Dispersion-curves for o=0.45. 


4, Conclusions ps ae 73 
(1) Ifc/v,>-1, the dispersion-curves always coincide with any one of the w= 1s 
curves given by (2.4), whatever value o may have. re aw 


(2) If c/v,>1, periods of the former and of the latter equation in (2. Ay 
should coincide with each other under the next condition, 


2m+ 1 = (21+ 1) (v,/v,) . 


ae Oo. Q2 03 O5 07 10 2 3 A oe de 
‘ager ta ===> 1 Vp/H 


Fig. 12. Dispersion-curves for o=0.35. ree 


ess Because odd orders of m must correspond to even orders of J, the lowest. 
order of m satisfying (4.1) is 1 for /=0, resulting in v,/v,=3.0 or o=0.44. 
If «<0.44, no curve indicated by m exists to the right of the curve 
indicated by /=0 with which M,®) always coincides in this case. 
ore If «>0.44, on the other hand, some curves indicated by m appear to the 
_ right of =0. The nearer o approaches to 0.50, the larger becomes the number 
of m-curves to the right of /=0. 
tal (3) On dotted lines in Fig. 10 to 14 amplitudes of displacement-potentials 
are zero, though the dispersion-curves may be calculated. 
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Fig. 13. Dispersion-curves for o=0.25. 


(4) It may be expected, on the other hand, that on chain-lines amplitudes 
will be large. Indeed amplitudes on full lines differ from zero, but the 
maximum amplitude must be attained on chain-lines at 


Tv,/H = 4/(21+1). (4.2) 
(5) It has been found to be wrong, if «>0.44 for M,,) and if «>0.22 for 
M,,"), to say that the lower the order of dispersive RAYLEIGH waves is, the 
larger may become the amplitude of the wave. This situation is completely 


different from the case of Love waves and of sound waves, consisting 
respectively of either S or P-waves alone. 


(6) (4.2) means “‘quarter wave-length law’. It was a wonder that such 


Fig. 14. Dispersion-curves for o=0.00. 


a law as (4.2) had been often experienced?) for dispersive RAYLEIGH waves in 
field-experiments, in spite of the existence of P-wave as well as of S-wave in 
the layer. But it has been understood in the present paper that among many 
branches of dispersive RAYLEIGH waves a few of them must correspond to 
cos s#=0 and others to cos BH=0. Amplitudes of the latter ones cannot be 
so large that the law containing v, alone will become important. It must be 
also noticed that the order of dispersive RAYLEIGH waves cannot be determined 
from the observed “wave-length law’. 

(7) Picking out M,"’, M,®) and the curve for /=0 alone from Figs. 10 to 
14, one will obtain Fig. 15. In this figure M,‘!) and M,") diverge from the 
curve for J=0 when o approaches to 0.50. Thus it should not be expected 
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Fig. 15. Comparison of dispersion-curves for the zeroth ra +i 
: order of M-waves to the curve for the zeroth Ra 


order of liquid waves. 


According to the present investigation, one sees, the process is very com- + 


plicated. When o approaches to 0.50, the zeroth order of liquid waves is "3 
constructed approximately of many higher orders M,,?)-waves and the first Ke 
order of liquid waves is approximately constructed of many higher crdens of : 


M,,1)-waves and so on. bat Pee 


- 5. Remarks | “ae 

If attention is confined in the region of Tv,/H<1 and c/v,~1 in Fig, 5, . 
a similar lattice to MINDLIN’s®) may be formed. The dispersion-curve in this 
region takes also “‘terrace-like structure’? named by MINDLIN in the case of ao 
plate. . 
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